JOURNAL OF GEOPHYSICAL RESEARCH, VOL. ???, XXXX, DOI:10.1029/, 



Modeling the Young Sun's Solar Wind and its Interaction with 
Earth's Paleomagnetosphere 

M. Glenn Sterenborg 

Department of Earth and Planetary Sciences, Harvard University, Cambridge, Massachusetts, USA 

0. Cohen 

Harvard Smithsonian Center for Astrophysics, Cambridge, Massachusetts, USA 

J.J. Drake 

Harvard Smithsonian Center for Astrophysics, Cambridge, Massachusetts, USA 

2 T.I. Gombosi 

Department of Atmospheric, Oceanic, and Space Sciences, University of Michigan, Ann Arbor, ML, USA 

Abstract. Wc present a focused parameter study of solar wind - magnetosphere inter- 
h— 5 action for the young Sun and Earth, ~3.5 Ga ago, that rehes on magnetohydrodynamic 

(MHD) simulations for both the solar wind and the magnetosphere. By simulating the 
quiescent young Sun and its wind we are able to propagate the MHD simulations up to 
Earth's magnetosphere and obtain a physically realistic solar forcing of it. We assess how 
<~~^ sensitive the young solar wind is to changes in the coronal base density, sunspot place- 

ment and magnetic field strength, dipole magnetic field strength and the Sun's rotation 
Z/y period. Prom this analysis we obtain a range of plausible solar wind conditions the pa- 

^ (-H leomagnetospherc may have been subject to. Scaling relationships from the literature sug- 

^ ^ gest that a young Sun would have had a mass flux different from the present Sun. We 

I evaluate how the mass flux changes with the aforementioned factors and determine the 

O importance of this and several other key solar and magnetospheric variables with respect 

> to their impact on the paleomagnetosphere. We vary the solar wind speed, density, in- 

^ terplanetary magnetic field strength and orientation as well as Earth's dipole magnetic 

field strength and tilt in a number of steady-state scenarios that are representative of 
young Sun-Earth interaction. This study is done as a first step of a more comprehen- 
sive effort towards understanding the implications of Sun-Earth interaction for plane- 
^ tary atmospheric evolution. 

o 

1. Introduction durations, up to several thousand years, the transitions be- 

tween polarities are thouglit to be accompanied by a strong 

The shape and size of Earth's magnetosphere is governed reduction of field strength, down to a tenth of its current 

by the interaction of the Earth's magnetic field and iono- strength. It is not yet fully understood how the field goes 

sphere with the solar wind. As such, the solar wind speed from one polarity to its opposite but paleomagnetic mea- 

^ and density, as well as the strength and orientation of both surements and dynamo models suggest that the tilt of the 

the Earth's magnetic field and the interplanetary magnetic dipole can experience excursions prior to reversal that can 

> field (IMF) play significant roles. Although the Earth's exceed 45° [Glatzmaicr and Roberts, 1995]. The field's mor- 

magnctic field at present has significant power in spheri- phology is also thought to change during such a transition, 

cal harmonic degrees beyond the dipole, wo may consider going ^om mainly dipolar to multipolar [Gubbins, 1994]. 

b it, on average, dipolar with its dipole axis coincident with All things being equal, such long -and short term changes 

^ Earth's rotation axis. However, the field may be considered i° magnetic field will have strongly affected the magne- 

vaxiable on many spatial and temporal scales. The Earth's , , , , . , , , 

J. , 4. 1, -4. -i J ij 1 the bun has also changed considerably over its liictime, 

dipole moment has seen its magnitude vary a great deal over r^ - i i ^^^r, ^ • , i n^r^m r. ■ , 

ii,! 4. rri J 4. 1 ronini u i, 4.1, 4. 4.1, Gudel, 2007: Guman and Englc, 2009 . It is assumed that 

Earth's history. Tarduno et al. 2010 have shown that Earth lu c v ^ • ^ j i ^ • i ^ ii, 

, J 4.- n u 1 n Al- /~i mi- • the Bun nas scon stages in its development similar to those 

had a magnetic neld as eaxly as ~3.45 Ga ago. iheir mea- u j j- ^ i ^ c 

, . ,. , , nil, ,1 ..^ rrnfT/ ,-.1 , observed tor younger stars analogous to a younger bun, on 
surements indicated a field strength 50-70% of the present ^^^^ main-sequence [see, e.g. Ribas ct al., 2005, 2010; Guinan 
field. This constrasts with a scahng relationship that sug- ^^^^ g^gj^^ 2009]. In keeping with main-sequence evolution, 
gests that the fie d strength scales with the rotation period ^^^^ g^j^j. luj^ynosity is assumed to have increased over its 
[Stevenson, 2003], which would suggest a stronger palco- ufetime, but the magnetic activity-related X-ray, EUV and 
magnetic field due to Earth's much shorter rotation period uV outputs have decreased owing to wind-driven magnetic 
in earlier times. Paleomagnetic measurements indicate the rotational braking and the associated decline in rotation- 
field has gone through many polarity reversals [McFadden powered dynamo activity [e.g. Skumanich, 1972; Guinan and 
and Merrill, 2000]. Irregularly spaced in time and of different Engle, 2009]. The solar wind is also highly variable, both on 

short -and long timescales. For a quiescent Sun we can make 
the distinction between the fast and slow solar wind in terms 
of low and high density, respectively. But for solar flares and 

Copyright 2013 by the American Geophysical Union. coronal mass ejections (CME), the solar wind speed and den- 

0148-0227/13/$9.00 sity can be much greater. Over longer timescales, current 
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observational evidence, albeit scant, indicates that main- 
sequence stellar wind mass flux declines with age [Wood 
et al., 2002, 2005]. This is important because stellar mass 
loss rates are directly linked to the magnetic braking that 
gradually slows stellar rotation and quenches dynamo and 
related magnetic activity. Another aspect of Sun-Earth in- 
teraction that is linked to the solar mass flux is the Faint 
Young Sun paradox, which could be resolved by the presence 
of greenhouse gases in Earth's atmosphere, a lower terres- 
trial albedo, or an increased mass loss rate of the young Sun 
[Sagan and Mullen, 1972; Guzik et al., 1987; Rosing et al., 
2010]. 

There is considerable topical interest in solar wind- 
magnetosphere (-ionosphere) interaction and in particular 
what this entails for Earth's magnetosphere during the early 
Archean (3.8 to 3.5 Ga), after the Late Heavy Bombard- 
ment and until the first appearance of fossil evidence for 
simple life [Kasting, 1993]. It is during this time that ex- 
ternal modification of Earth's atmosphere (e.g. by impacts) 
became infrequent enough to allow for a long-term steady 
state in the physical and chemical properties of the atmo- 
sphere. Furthermore, this is as far back in Earth's history 
that we have evidence of the presence of a planetary mag- 
netic field [Tarduno ct al., 2010]. 

Understanding the magnetospheric response to variations in 
solar output is a necessary first step in understanding how 
these variations ultimately affect Earth's atmosphere and its 
evolution. This field is also of interest in the larger context 
of interaction between (magnetized) extrasolar planets and 
their parent stars and its implications for the development 
of planetary atmospheres. 

In this study we focus on how the young Sun interacted 
with Earth's paleomagnetosphorc. Characterization of the 
paleomagnetospliere, in terms of parameters such as rnag- 
iictopausc stand-off distance, flank distance, polar cap size 
etc. is often carried out by scaling relationships [Siscoe and 
Chen, 1975; Saito et al., 1978; Vogt and Glassmeicr, 2001]. 
More recently, nurncrial MHD simulations have been em- 
ployed to build on this work [Ziegcr ct al., 2006a], in which 
several scaling relationships were investigated between the 
Earth's paleomagnetosphere, its dipole moment and the in- 
terplanetary magnetic field. Zieger et al. [2004] studied the 
case of an equatorial dipolar paleomagnetosphere in a large 
parameter space study by varying the strength of Earth's 
dipole moment, its dipole tilt and the IMF also using MHD 
simulations. A similar study looked at quadrupolar pale- 
omagnetospheres as representative of rrmltipolar fields dur- 
ing geomagnetic polarity reversals [Vogt et al., 2004]. Zieger 
et al. [2006b] researched magnetosphere-ionosphere coupling 
for the paleomagnetosphere in terms of the transpolar po- 
tential and associated field aligned currents. It is noteworthy 
that many of these studies were carried out for constant solar 
wind speed and density. In related work but focusing on hot 
Jupiter-sized exoplanets, Griefimeier et al. [2004] used rela- 
tively simple scalings for the planetary dipole moment, core 
radius, magnetopause stand-off distance and stellar wind as 
well as a parametric model of the magnetosphere, to study 
their magnetospheric and atmospheric evolution. 

Rather than investigating a large parameter space that 
is expected to include most Sun-Earth interaction scenar- 
ios, from weak multipolar magnetic fields during polarity 
reversals to strong steady axial dipoles, we attempt a more 
focused approach aimed at simulating most likely scenarios 
of young Sun-Earth interaction. We start by numerically 
modeling the young Sun and obtain a range of likely values 
for the solar wind speed and density. These then serve as 
a forcing of an MHD model of Earth's paleomagnetosphere. 
Magnetospheric response to solar output is cast in terms of 



several key observables such as the subsolar magnetopause 
stand-off distance, the magnetopause flank distances, polar 
cap extent and the magnetic and plasma pressures inside 
the magnetosphere. 

The main contribution of this study to the field of solar 
wind-magnetosphere interaction comes from our characteri- 
zation of the young Sun and obtaining physically jusitifiable 
solar wind scenarios with which to force the Earth's magne- 
tosphere. This paper, which is readily subdivided in several 
parts, is structured as follows: In section 2 we discuss our 
approach to modeling Sun-Earth interaction using the Space 
Weather Modeling Framework, in particular how we model 
the young Sun and Earth's magnetosphere. We carry out 
a sensitivity analysis to determine how responsive the solar 
wind is to changes in sunspot placement, sunspot magnetic 
field strength, dipole magnetic field strength and rotation 
period in section 3, and in the process obtain solar wind so- 
lutions most representative of the young Sun. In section 4 
we assess how the solar mass flux from our models compare 
with literature estimates. The MHD simulations of Earth's 
magnetosphere forced by these solar wind solutions are dis- 
cussed in section 5. We close with our conclusions in section 
6. 

2. Approach 

For our simulation of the young Sun and its influence 
on Earth's paleomagnetosphere we use several modules of 
the Space Weather Modeling Framework (SWMF). This 
is a high-performance flexible framework for physics-based 
space weather simulations, as well as for various other space 
physics applications. The SWMF is able to integrate numer- 
ical MHD models of the solar corona, solar eruptive events, 
the inner heliosphere, solar energetic particles, the global 
magnetosphere, the inner magnetosphere, the Earth's radi- 
ation belt, ionosphere electrodynamics, and the upper at- 
mosphere into a single coupled system. The SWMF enables 
comprehensive simulations that are not possible with the in- 
dividual physics models. For a full description of the frame- 
work we refer the reader to [Powell et al., 1999] and [Toth et 
al, 2005]. The SWMF delivers the capability to assess the 
entire system from Sun to Earth in a self-consistent man- 
ner and allows us to solve for the MHD variables in Earth's 
magnetosphere having started from the solar corona. 

Our approach is straightforward in that we start at the 
Sun and work our way towards the Earth's magnetosphere. 
We employ observations of the Sun to constrain a solar 
corona model and modify it to represent the young Sun. We 
use this model's output to force the Earth's magnetosphere, 
the model of which is also modified to represent a paleo- 
magnetosphere, and obtain a physically realistic response. 
To account for the myriad uncertainties that are inherent in 
this problem, we vary our input parameters over a range of 
realistic, physics based, values. We assess the response by 
inspecting such outputs as subsolar magnetopause stand- 
off distance Rmp, magnetopause flank distances Rdd (dawn- 
dusk) and R„B (north-south), polar cap latitudinal extent 

Opc Hilld 3>rG3j .ApQ. 

3. Simulating the young Sun 
3.1. Theory 

We use the solar corona (SC) model of the SWMF [Cohen 
et al., 2007, 2008] to simulate the (young) Sun. This model 
is based on the generic MHD BATS-R-US model [Powell et 
al., 1999] and is driven by surface magnetic field maps which 
determine the initial potential magnetic field and are used 
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to scale the boundary conditions on the solar photosphere. 
The energy deposition to the solar wind is based on the ob- 
served inverse relation between the amount of expansion of a 
particular magnetic flux tube and the terminal speed of the 
wind flowing along this tube. For an initial potential field 
distribution (which we calculate using harmonic coefficients 
supplied by the Wilcox Solar Observatory which observes 
the Sun's photospheric magnetic field), Wang and Sheeley 
[1990] have determined a flux tube expansion factor as: 



Rss 



Bq{Rss 



(1) 



where Rq and B{Rq) are the solar radius and the mag- 
nitude of the magnetic field at the flux tube base surface, 
respectively. B{Rss) is the magnetic field magnitude of the 
same flux tube at a height Rss, which is the height of the so- 
called "source surface". The potential field approximation 
assumes by definition that at the source surface, all field lines 
are fully radial and open; the source surface is commonly set 
at Rss = 2.5Rq. Based on the above definition of fs, Arge 
and Pizzo [2000] have derived an empirical relation between 
the expansion factor and the observed solar wind speed at 
~1AU (essentially the terminal speed) u,,„,(r — > oo,l//,). 
This empirical model is known as the Wang-Shcclcy-Argc 
(WSA) model. 

In the MHD model used here, the energy necessary to 
power the solar wind is determined using the empirical re- 
lation above in the following way. Assuming that Usw is 
known, we can also assume that the total energy of the so- 
lar wind far from the Sun equals the bulk kinetic energy of 
the plasma, Us^/2. On the solar surface, the total energy 
equals the enthalpy of the plasma minus its gravitational 
energy. Assuming the conservation of energy along a partic- 
ular fiux tube, we can relate the two ends using Bernoulli's 
integral, where we have assumed an adiabatic expansion and 
a zero speed on the solar surface: 



70 kbTo GMe 



70 — 1 mp 



R, 



(2) 



Here kb is the Boltzmann constant, rnp is the proton mass, 
G is the gravitational constant. To is the boundary condi- 
tion value for the temperature, and 70 is the surface value of 

the ratio of specific heats. Equation 2 enables us to specify 
tfie value of 7 at the coronal base as a function of the final 
speed of the wind that originates from this point. The value 
of 7 is close to unity near the solar surface but is closer to 
1.5 at 1 AU [Totten et al., 1995]. A gas with 7^1 can 
be considered to have a higher internal (turbulent) energy 
than a gas with 7 ~ 1.5. Therefore, this observed change in 
the value of 7 can be associated with the amount of inter- 
nal energy that is released to accelerate the solar wind [see 
Roussev et al., 2003, for a full description]. Based on the 
surface value of 7 (determined by Equation 2), we can define 
a volumetric energy source term, E~,.{r,jo), which vanishes 
as 7 — >■ 1.5. With E^y and setting an inner boundary condi- 
tion for the density we now self-consistently solve the set of 
MHD equations 
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+ V • (uB - Bu) = 



which couples the conservation of mass, momentum, and 
energy laws as well as the magnetic induction law where F 
represents a constant 7 = 1.5. The V • B = condition is 



enforced using the 'eight- wave' method [Powell et al., 1999]. 
We thus obtain a steady-state solution for the solar wind 
speed, density, pressure and magnetic field strength. In do- 
ing so we generate a solar wind solution that is consistent 
with the distribution of the Sun's surface magnetic field and 
is constrained by the empirical WSA solar wind model. 

In all the SC simulations presented here, we use a non- 
uniform grid in the SC model (0 — 24i?Q) with grid size of 
Aa:: ~ W^^Rq near the solar surface. The grid is also dy- 
namically refined near the coronal and hcliospheric current 
sheet in order to improve the solution. In the Inner Helio- 
sphere (IH) module {17 Rq — lAU) the grid size near Earth 
is Aa; m 0.5Rq. SC simulations were performed using the 
PLEIADES system at the NASA Advanced Supercomput- 
ing (NAS) Division. 

3.2. Implementation 

For the purposes of modeling the solar wind and corona, 
the most relevant feature of a young, 'active' Sun is the pres- 
ence of a stronger surface magnetic field than wc observe to- 
day. Simple considerations predict that the efficiency of the 
magnetic dynamo scales with tQ, where r is the convective 
turnover time in the layers where the dynamo is operating, 
and n is the angular velocity of these layers [e.g. Durney and 
Latour, 1978]. The expected trend of increasing magnetic 
field strength with increasing rotation rate is borne out, at 
least qualitatively, by magnetic field observations of solar- 
like stars [e.g. Petit ct al., 2008], and the behaviour of mag- 
netic proxy indicators such as X-ray, EUV and UV emissions 
[e.g. Guinan and Engle, 2009; Pevtsov et al., 2003]. Petit 
et al. [2008] have used the Zeeman Doppler Imaging tech- 
nique to map out the surface magnetic field of Sun-like stars 
with different rotation velocities. They find an increase in 
the mean stellar magnetic field with rotation velocity, from 
3.6±1 G in HD 146233 with a rotation period of 22.7 days to 
42±7 G in HD 73350 with a period of 12.3 days. HD 146233 
is then not only very similar to the present-day Sun in ro- 
tation period, but also in terms of its global magnetic field. 
At longer rotation periods, poloidal fields (mainly dipolar) 
completely dominate the global field, as we see on the Sun. 
However, Petit et al. [2008] find that toroidal components 
increase with decreasing period such that 50% of the mag- 
netic flux of HD 73350 resides in toroidal flelds. We use these 
magnetic field measurements to guide our simulations. In 
order to mimic the surface magnetic field of the young Sun, 
wc take detailed observations of the present-day Sun and 
modify the magnetic field strength and surface distribution 
of these in ways indicated by observations of young solar 
analogs. 

The baseline for each simulation is a high reso- 
lution Michelson Doppler Imager (MDI) magnetogram 
(http://sun.stanford.edu/synop/) for Carrington Rotation 
(CR) 1958 (taken in January 2000 during solar maxiirmm). 
During CR1958 the Sun was relatively active with a large 
number of active regions on the disk, unlike solar minimum 
conditions during which the disk is clear of sunspots and 
the solar magnetic field is nearly dipolar. We 'modify' so- 
lar activity in two ways. First, we decompose the surface 
magnetic field in terms of its weak component, which repre- 
sents its dipolar component, and a strong field component, 
i.e. the sunspots' magnetic field. The threshold value for 
this decomposition is 10 G. This allows us to investigate the 
role each plays by independently scaling them. Secondly, 
by moving sunspots to higher latitudes we obtain artifical 
magnetograms with a realistic spot distribution as observed 
for young stars such as AB Doradus [Hussain et al., 2002). 
Given that the young Sun is assumed to have a faster ro- 
tation rate than the present Sun [Giidel, 2007; Guinan and 
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Engle, 2009], we carry out simulations with the present ro- 
tation period of 27 days and with a shorter period of 15 
days. 

3.3. Sunspot latitude sensitivity analysis 

Before we attempt to simulate the young Sun solar wind 
we first ascertain how well the SC model performs for 
the case of the present Sun. We extract the SC solution 
along the trajectory of the Advanced Composition Explorer 
(ACE) spacecraft, which is situated at the LI Earth-Sun La- 
grange point, at ~1 AU, and compare it with in-situ ACE 
data taken from http://cdawob.gsfc.nasa.gov. Wo repeat 
this procedure for each simulation and obtain solar wind so- 
lutions over the course of a full solar rotation. Simulations 
in subsequent sections follow the same method. 
We also test sensitivity of the solar wind solutions, i.e. of so- 
lar wind speed, density, magnetic field at 1 AU, to latitudinal 
placement of sunspots in the modified magnetogram. Figure 
1 shows the original and modified magnetograms used in this 
sensitivity analysis. We carry out several simulations with 
four magnetograms: 'shiftOO', which represents the original 
unmodified magnetogram for CR1958; 'sliiftSO', where the 
sunspots have been moved to ~30° latitude; 'shiftGO', with 
the sunspots at ~60° latitude and 'CR2074-shift00', which 
represents a comparison magnetogram for a solar minimum. 
For these cases, both sunspot and dipole magnetic fields are 
scaled by a factor of two, which is a baseline scaling always 
applied to reduce discrepancies between total magnetic flux 
predicted by the magnetogram driven potential field model 
and that observed at 1 AU [Cohen et al., 2008]. Comparing 
the simulated solar wind to ACE data for the present Sun, 
Figure 2 shows the predicted and observed solar wind con- 
ditions: plasma speed u (km/s), magnetic field strength B 
(nT), plasma immbcr density n (cm^'^) and the local value 
of the mass flux pu {Mq yr^^ m^"^) at 1 hour averages along 
the location of Earth during CR1958. 

As is clear from the ACE data the solar wind is quite vari- 
able. The record reveals two corotating interacting regions 
(CIR) as well as a coronal mass ejection (CME). CIRs are 
compression regions created when the fast solar wind ex- 
ceeds the slow solar wind as the two propagate into inter- 
planetary space. In general, two CIRs appear in a single 
CR, and they arc associated with a smooth increase in wind 
speed (transition from slow to fast wind) , and an increase in 
B, n, and T due to the compression in the CIR. A CME is 
a solar eruption that carries with it large amount of ionized 
gas and is propagating through interplanetary space with a 
shock at its front. Therefore, a signature of a CME in the 
in-situ measurement will appear as a sharp increase in the 
plasma parameters, see the end of the record. There are 
some other indicators, such as composition data, for a CME 
in the ACE data that are not discussed, since here we are 
interested in the steady-state solution for the solar wind and 
modeling of CMEs is omitted. 

For each observable we also show the average value for the 
entire year 2000. In the understanding that we are compar- 
ing a steady-state solution for the quiescent Sun to a single 
realization of the Sun's output we note that the solar wind 

solutions for the original magnetogram, 'shiftOO', perform 
very well in capturing the main features with the exception 
of the CME which we do not model. 

In order to safeguard against a particularly (un)fortunate 
comparison of the solar wind solution to the ACE data we 
also show how the solution varies with latitude at 1 AU. 
These are the 'shiftOO-up' and 'shiftOO-down' solutions which 
were taken at ±15° away from the ecliptic plane. The so- 
lution is well behaved and does not vary significantly away 
from Earth's orbit. For each panel the solar wind solutions 
appear to slightly lead in phase. This is likely due to a slight 



overcstimation of speed causing features to arrive at 1 AU 
a little sooner. For each variable the solar wind solutions 
also perform very well in terms of magnitude. We note that 
we slightly underestimate the magnetic field at 1 AU, how- 
ever this is due to how we model the solar corona. Since 
we start from a specification of the magnetic potential field, 
any observed open magnetic flux originating from outside 
the polar coronal holes, e.g. at lower latitudes and around 
active regions, goes unmodelled yielding a discrepancy in 
magnetic field [see Cohen et al., 2008]. 

The mass flux panel clearly shows that density is the dom- 
inant factor with the velocity only acting as a prefactor. 
As we move the sunspots to increasingly higher latitudes, 
'shiftSO' and 'shift60', we sec that the speed and density 
are quite sensitive to latitudinal sunspot placement. Fig- 
ure 2 and the magnetic field less so. For a 30 degrees shift 
the features in each panel are less pronounced and for a 60 
degrees shift they have all but disappeared, retaining only 
the average magnitudes. Figure 3 and Table 1 show how 
these solutions vary in terms of mass flux. We see that with 
increasing latitude of the sunspots the mass flux increases 
slightly but the distribution remains roughly unaffected. 
For the sake of completeness we also show a solution, 
'cr2074-shift00', representing a solar minimum. This solu- 
tion seems able to capture some of the larger scale features 
in the data but, as it should do, under- and overestimates 
the solar wind speed and density respectively. 
What these simulations show is that the SC model of the 
SWMF is very capable of modeling the solar wind for the 
present Sun, which makes us confident that it will be able 
to reasonably model the young solar wind as well. 

We also carried out simulations for the magnetogram ob- 
tained during CR2010 and obtained consistent results, but 
with slightly poorer agreement with observations. This was 
due to an uncharacteristically active Sun for that solar max- 
imum, which featured the 2003 Halloween storms, and the 
associated "contamination" of the more steady solar wind 
components with CMEs that are not included in the model. 

3.4. The role of coronal base density 

The solar wind parameters derived from the MHD mod- 
els presented here are sensitive to the assumed value of the 

coronal base gas density. The density of gas in the different 
regions of the present-day solar corona is well-established 
from extreme ultraviolet spectroscopy, and ranges from 10*^ 
cm^'^ in the quiet Sun and diffuse limb to typically ~5 X lO" 
cm^'^ for active regions [e.g. sec Keenan et al., 2008; Shcstov 
et al., 2009; Keenan et al., 2010, for recent assessments]. 
Spectra of the full-disk Sun, like spatially-unresolved stellar 
observations, largely refiect the properties of the brighter, 
more dense active regions, and Laming et al. [1995] found 
densities of 1-3 x 10^ cm~^ from the full-Sun EUV spectrum 
of Malinovsky and Heroux [1973]. 

The behaviour of plasma density in solar-like stars indicates 
that the more active Gyr-old Sun is likely to be character- 
ized by somewhat higher densities. X-ray spectra from the 
Chandra and XMM-Newton satellites show plasma densities 
in solar-like activity stars, such as Procyon and a Cen A 
and B, are similar to full-disk solar measurements at a few 
lO'' cm"'^ [e.g. Mcwe et al, 1995; Schmitt ct al., 1996; Drake 
et al., 1997; Ness et al., 2002]. Towards higher activity level 
stars, densities transition to a few 10^° cm~^ at base coro- 
nal temperatures of 1-2 x 10*^ K [e.g. Testa et al., 2004; 
Ness et al., 2004]. These measurements are based on X-ray 
emission from mostly closed magnetic loops, and it is not 
yet certain to what extent these properties and trends are 
shared by open field regions from which the wind largely 
emanates. Nevertheless, the indication is that the coronal 
base density for the ~ 1 Gyr Sun would have been about 
a factor of ten larger than it is today. Figure 2 shows the 
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solar wind solution for a model in which we specify a high 
coronal base density of 10^" cm~^ as compared to the oth- 
erwise used 5 X 10* cm~^: 'shiftOO-hd'. This solution for a 
zero degree shift of the sunspots shows a much smaller so- 
lar wind speed and concomitantly much higher density. As 
such it reveals how the solar wind responds to an increase in 
coronal base density and yields a possible young Sun-Earth 
scenario. Note the large increase in mass flux for this solu- 
tion, Table 1. 

Sensitivity of the SC model to longitudinal sunspot place- 
ment, i.e. to choice of magnctogram, was addressed in Co- 
hen ot al. [2008] where SC output was compared to one year 
of ACE data (CR1972-CR1984) during which period the Sun 
was experiencing solar maximum conditions. The model was 
shown to perform well and captured the salient features of 
the solar wind such as the periodic corotating interacting 
regions and transients from slow to fast solar wind. As such 
we are confident our choice of magnetogram does not bias 
our solar wind solutions. 

3.5. Sunspot and dipole magnetic field 

Let us now look at solar wind solutions for which we in- 
dependently vary the sunspot magnetic field, 'sf, and the 
background dipole magnetic field, 'df '. In general, the dipole 
component of the solar magnetic field generates fast wind 
originating from regions of open field lines at high latitudes 
(coronal holes), and slow wind originating from open field 
lines at the boundary of the large helmet streamers. On 
top of this ambient, large-scale, bi-modal wind structure, 
smaller scale modulations appear due to solar wind that 
comes from open field lines at the vicinity of the spots. The 
mass flux from these regions is usually higher than the mass 
flux associated with the high-latitude fast wind since the 
source regions near the spots arc denser than the ambient 
coronal holes. In the case of surface field dominated by the 
dipole component, the mass flux distribution is dominated 
by the bi-modal ambient wind structure so that the mass 
flux distribution is bi-modal as well. In the case of strong 
spots however, local variations in the mass flux distribution 
due to the spots arc bigger. 

Figure 4 shows the magnetograms that were used as inputs 
to the SC model. For each magnctogram wc placed the 
sunspots at 60° latitude and scaled the sunspot and dipole 
magnetic fields. We use the nomenclature 'sfxdfy' where 'x' 
and 'y' are scaling factors. In Figure 5 we show solar wind 
solutions for these magnetograms. All solutions are for a 
rotation period of 27 days, with the exception of those that 
have '15d' in their name, which are for 15 days. Again we 
plot u, n, B and mass flux and show ACE data for reference. 
For the first three solutions, 'sf2df2','sf2df5' and 'sf2dfl0', 
we only increase the dipole magnetic field which, on aver- 
age, results in an increase in speed and a slight increase 
in plasma density. These scaling factors were guided by 
the Petit et al. [2008] magnetic field maps for solar analogs 
mentioned above, that indicate a higher poloidal field by an 
order of magnitude than the current Sun for a shorter rota- 
tion period of 12 days that is more appropriate to the young 
Sun modeled here. The magnetic field strength stays ap- 
proximately the same as does the local mass flux. Looking 
at the mass fluxes for those runs, see Figure 6 (left column) 
and Table 1, we note that the mass flux slightly decreases 
with increasing dipole strength and has a tendency to de- 
crease near high latitudes and concentrate at lower latitudes 
at 1 AU. 

Next we vary only the sunspots, increasing the spot mag- 
netic field by factors of 2, 5 and 10 (models 'sf2df2', 'sf5df2' 
and 'sfl0df2'; see Figure 4). The contours in this flgure are 



saturated to clarify the increase in magnetic field. The mag- 
netic field strength for 'sf5df2' and 'sfl0df2' reaches several 
kC, in accordance with observations of stellar magnetic field 
for young solar analogs [Giidel, 2007]. As Figure 6 (right 
column) and Table 1 indicate, at 1 AU this increases the so- 
lar wind speed more severely, the magnetic field intensifies, 
the density slightly decreases and the local mass flux sees 
slightly stronger variations. 

The mass flux is enhanced by increasing the sunspot mag- 
netic field and becomes concentrated at larger local regions 
comparing to the strong dipole cases. The wind speed in- 
creases for 'sf5df2' but decreases for 'sfl0df2'. This is due 
to the fact that when the high-latitude spots become strong 
enough, they close down high-latitude open field lines. This 
process eliminates much of the fast wind originated from 
coronal holes and associated with lower mass fiux, and fill 
the interplanetary space with wind associated with the vicin- 
ity of the active regions. 

As noted earlier, because solar-like stars lose angular mo- 
mentum over time to their winds, they will have had shorter 
rotation periods in the past. We tested the impact of this by 
repeating the last solar wind model with a 15 day rotation 

period, 'sfl0df2-15d', the main result of which is that the 
solution is slightly phase shifted with respect to the same 
solution for 27 days, 'sfl0df2'. 

The solar wind solutions 'sf2dfl0_15d_hd' and 'sfl0df2_15d_hd' 
feature a strong dipole or sunspot magnetic field respec- 
tively, a short rotation period and a high coronal base den- 
sity in an effort to determine their combined effects on the 
solar wind and obtain a realistic possible scenario for the 
young Sun. In keeping with the solution 'shiftOO-hd' of the 
previous section. Figure 5 and Table 1 show a relatively 
slow solar wind, high wind density, a stronger magnetic field 
and naturally a high (local) mass flux compared to similar 
non-hd solutions. Figure 6 also shows this on the bottom 
row where we point out the higher mass flux limits on the 
colorbar. Although the mass fiux distribution for the 'hd' 
solutions does not reveal any noteworthy specific features, 
it is clear that the coronal base density is a dominant term 
controlling the solar wind mass flux. 

We note that all solutions shows a strong anti-correlation be- 
tween speed and density as expected [McComas et al., 2007). 

In summary, we have parameterized the SC model to ob- 
tain a reasonable agreement with present 1 AU data. Obser- 
vations of young stars reveal that stellar spots are located 
at higher latitudes [Hussain et al., 2002] and so we have 
shifted the spots by 30° and 60° to determine the effect on 
the mass flux distribution. We also used the SC model to 
test how variations of the magnitude of the sunspots' mag- 
netic field, as well as the ambient weak field affect the mass 
fiux distribution. In so doing wo have covered a compre- 
hensive set of scenarios as far as the surface magnetic field 
distribution is concerned. This also gives us constraints on 
the mass fiux, which is based on a detailed model parame- 
terization instead of on scaling relationships such as those 
by Wood et al. [2002, 2005]. As mentioned before, Ziegcr 
et al. [2004]; Vogt et al. [2004]; Zieger et al. [2006a] carried 
out comprehensive parameter studies in which they mainly 
varied dipole moment magnitude and tilt but, with the ex- 
ception of the IMF, kept the solar wind parameters constant. 
Here we have obtained a range of solar wind scenarios which 
we can now use to force Earth's paleomagnetosphere with. 

4. Other mass flux indications 

Wood et al. [2002, 2005] used Lya spectra to infer hot 
hydrogen absorption in astrospheres - the stellar analogs of 
the heliosphere - which they used diagnostically to determine 
stellar mass loss rates. To quantitatively measure mass-loss 
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rates based on astrospheric absorption requires the use of 
hydrodynamic models of the astrosphere, which in turn are 
based on a model of the Sun's heliosplicre that reproduces 
licliospheric HI absorption. This requires assumptions re- 
garding: solar wind velocity, proton density and tempera- 
ture, in addition to local interstellar medium (ISM) prop- 
erties. Application to other stars relics on ISM properties 
not varying significantly between the different stars studied. 
Furthermore, an additional, and perhaps dominant source of 
error in HI absorption measurements comes from the need 
to assume a form for the unabsorbed stellar Lya profile. 
Wood et al. [2002, 2005] estimated a relationship between 
the stellar mass loss rate M and X-ray flux Fx 



M oc Fi34±o.i8 



(4) 



This can be combined with estimates from Ayres [1997] of 
rotation rates, Vrot, for stellar analogs [Wood et al., 2005] 



Vrot oc t 



(5) 



which shows that as stars age their rotation rates become 
smaller, due to magnetic braking, and 



Fx oc V„, 



(6) 



which illustrates that a smaller rotation rate means less mag- 
netic activity and concomitantly less X-ray flux. The com- 
bination of equations 4, 5 and 6 yields 



M oc f 



-2.33±0.55 



(7) 



Using equations 4 and 7 Tarduno et al. [2010] arrived at a 
range of mass loss rates between 2.4xl0~^^ to 1.5xl0~^^ 
Mq yr~^ for the Archean Earth. 

Holzwarth and Jardine [2007] used a different approach, 
based on a magnetized solar wind model assuming power 
laws for the dependence of thermal and magnetic wind pa- 
rameters on the stellar rotation rate, to arrive at a mass flux 
for cool main sequence stars of no more than 2xl0~^^ Mq 
yr~^, i.e. less than 10 times the present solar mass loss rate 
of 2.1x10"" Mq yr-i [Feldman et al., 1977; Giidel, 2007]. 
With the exception of the solar wind solutions with high 
coronal base densities, 'shift00_hd', 'sf2dfl0_15d_lid' and 
'sfl0df2_15d_hd', which have mass fluxes between ~ 1 — 1.5 x 
10~^^ Mq yr^"', our solutions have mass fluxes within the 
range of 5 - 8x10^""^^ Mq yr^^ , well below the literature 
estimates by an order of magnitude. From this we may in- 
fer that cither scaling relationships such as those derived 
by Wood ct al. [2002, 2005] and Holzwarth and Jardine 
[2007] are overestimating young solar mass flux, or that 
higher coronal base densities are warranted in our models 
of the young Sun's solar wind. Given the expectation that 
younger solar-like stars are characterized by higher coro- 
nal base densities, see subsection 3.4, the latter seems more 
likely. However, we also note that the astrospheric response 
to stellar winds on which the Wood et al. [2002, 2005] mass 
loss estimates are based has a timescales of years or more. 
As such, these estimates include the average effect of mass 
loss through CMEs. The partition of mass loss between the 
steady wind and more impulsive CMEs is not known for 
stars significantly more active than the Sun, but based on 
the higher flaring rate of young solar analogs it is reasonable 
to expect that a greater proportion of the mass loss could 
occur impulsively. It is plausible that the difference in mass 
loss rates for our steady wind models and the estimates of 
Wood et al. [2002, 2005] is due to the impulsive component. 

5. Forcing Earth's paleomagnetosphere 

In order to assess the effect of various solar wind scenar- 
ios on the Earth's paleomagnetosphere we rely on the Global 



Magnetosphere (GM) MHD model of the SWMF, the input 
parameters of which we modify to siirmlate a paleomagneto- 
sphere. Based on the solar wind solutions of section 3, using 
the time-averaged values of Table 1, we develop several sce- 
narios with which to simulate Earth's paleomagnetospheric 
response to the young Sun's possible forcing (see Table 2). 
The GM domain ranges from about 30 Re on the daysidc 
to about 220 Re on the nightside and about 130 Re per- 
pendicularly to the Sun-Earth line. The inner boundary of 
the GM is 2.5 Re from the center of Earth with its bound- 
ary conditions set by the Ionospheric Electrodynamics (IE) 
model of the SWMF, [Ridley and Liemohn, 2002; Ridley, 
2005]. We use three levels of grid reflnement in the GM 
model with a minimum grid size of 1/8 Re closest to the 
Earth and a maximum grid size of 8 Re- 

Although the Earth's magnetic field is not strictly dipo- 
lar, over geologic timescales it does average to a dipole field 
with its axis of symmetry coincident with Earth's rotation 
axis [Merrill and McElhinny, 1983]. As the magnetic field 
reverses polarity its morphology changes as does its dipole 
axis which can experience excursions exceeding 45° . In this 
study we assume a dipolar magnetic field and a number of 
dipole tilt scenarios. 

Tarduno et al. [2010] have made paleomagnotic measure- 
ments of 3.5 Ga old samples that indicate a virtual dipole 
moment (VDM) ~50-70% of the present dipole moment, 
but representing field averages spanning only decades to cen- 
turies. We must be cautious when using VDM's as represen- 
tative of actual dipole moments, especially when the mea- 
surements the VDM is based on all come from the same, or 
similarly placed, sites, in which case the spatial distribution 
of the measurements is not adequate to separate the dipole 
from non-dipolar components. As such, VDM's may under- 
or overestimate the actual dipole moment. 
Given the geomagnetic reversal rate over the last 160 Ma of 
~5 per million years, with an average reversal duration of 
10^ - 10" years [McFadden and Merrill, 2000] and assuming 
this reversal rate has not changed significantly over the last 
3.5 Ga we estimate that the (purely statistical) probability 
of the aforementioned measurements of sampling normal, 
non-transient, field conditions is between 95-99.5%. How- 
ever dynamo theory suggests a scaling relationship that in- 
dicates that magnetic field strength inversely scales with the 
square root of the rotation period, implying a stronger pa- 
leomagnetic field. 

Planetary dynamos are believed to exist in a regime where 
the Coriolis and Lorentz forces are balanced by each other, 
which is captured by the Elsasser number 



(8) 



a ratio of the two, and should always be of order unity 
[Stevenson, 2003]. Here a is the electrical conductivity and 
p the density. We can obtain a scaling relationship between 
the magnetic field strength B and rotation period T = 27r/f2, 
where f2 is the rotation rate. 



SocT 



-1/2 



(9) 



Hansen [1982] suggests a rotation period of ~ 15 hours 
for the Archean Earth, due to the closer proximity of the 
Moon, which implies a ~30% increase in the magnetic field 
strength. Here we account for both possible scenarios: a 
stronger dipole field that scales with the rotation period 
and a weaker field as indicated by paleomagnetic measure- 
ments. 
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Wc assess each GM solution in terms of several observ- 
ables. The subsolar rnagnetopause stand-oflF distance, Rmp, 
is determined by the position of mininmm plasma speed 
at the stagnation point on the Sun-Earth line. The magrie- 
topause dawn-dusk flank distance, Rdd, and the north-south 
flank distance, Rns , are determined by finding the maximum 
current density on the magnetopause. We assess the extent 
of the polar cap by inspecting the boundary between open 
and closed field lines at the inner boundary of the simulation, 
the IE grid. We find the polar cap magnetic co-latitude, Opc, 
for the dayside and nightsidc. We also calculate the polar 
cap area relative to that of a GM solution for the present 
Earth. 

In scenarios A-C in Table 2 we vary Earth's dipole mag- 
netic field from 25-133% of the present field. Solar wind 
speed and density are based on the young Sun solutions 
with strong sunspot magnetic fields, 'sfl0df2' and 'sfl0df2- 
15d'. Figure 7, top-row, shows how the paleomagnetosphere 
changes with decreasing terrestrial dipole field strength. 
Qualitatively, it is clear that the magnetopause stand-off 
distance strongly decreases and the magnetosphere becomes 
very compressed. Quantitatively, Table 2 shows that Rmp 
more than halves, with strongly decreasing Rdd and Rna and 
increasing polar cap. 

During a geomagnetic polarity reversal the dipole field 
can exhibit several transitions, [Zieger et al., 2004; Gubbins, 
1994]. The field may become strongly multipolar before it 
settles into a dipole field of the opposite polarity, the dipole 
field may severely weaken and gain in strength with the op- 
posite polarity or the dipole field may experience several 
excursions, where the dipole axis tilts, before it finally flips. 
We address this latter possibility in scenarios D-F, which 
are based on the first four solutions of Table 1, 'sf2df2' - 
'sf5df2'. The middle-row of Figure 7 shows the response to 
these scenarios and from Table 2 it is clear that while Rmp, 
Rdd and Rna do not change much with dipole tilt, although 
-Rns does become asymmetric, the polar cap becomes much 
smaller. 

In scenario G we run a similar scenario to D-F but for 
a magnetic fleld similar to that measured in Tarduno et al. 
[2010]. We obtain a Rmp of 6 Re which is just within the 
upper margin of errror of their estimate. 

Earth's magnetosphere is not constantly exposed to a 
southward IMF along with the subsequent magnetic recon- 
nection and magnetic field erosion. The effect of a northward 
IMF is shown with scenario H where we see that the mag- 
netopause stand-off distance is relatively greater, 6.25 Re, 
than that for the previous scenario G, but the main effect 
is that there is no polar cap of any significance, i.e. the ab- 
sence of reconnection prohibits the solar wind from directly 
accessing the upper atmosphere. We have also run a north- 
ward IMF scenario for a strong dipole field, not shown here, 
which yielded Rmp '-^9 Re- 

Scenarios I through M, based on solar wind solutions 
'shiftOO-hd' (I), 'sf2dfl0-15d-hd' (J, K) and 'sfl0df2-15d-lid' 
(L, M), investigate the effects of a higher coronal base den- 
sity, which is characterized by a slow solar wind, a high 
wind density and a large mass flux (compared with present 
day). We also vary Earth's magnetic fleld between 50 and 
133 % of the present fleld. The effects of these scenarios on 
the magnetosphere are stand-off distances between 4 and 6 
Re depending on Earth's magnetic fleld strength, flank dis- 
tances between 5 and 9 Re, very high plasma pressures and 
relatively larger polar caps, see Figure 7 which illustrates 
this with scenarios J, K and L, as well as Table 2. However, 
this ~ 15 — 50% increase in polar cap is much less compared 



with polar cap estimates of Tarduno et al. [2010]. They used 
a scaling relationship from Siscoe and Chen [1975] 

cos(V)=(^J P^^''cos(e;,) (10) 

where is the present latitude of the division, P is the solar 
wind dynamic pressure normalized at the magnetopause to 
its present value of ~2 nPa and M, M' are the Earth's dipole 
moments in the Archean and the present respectively. Tar- 
duno et al. [2010] suggested a polar cap increase of ~300% 
with respect to present conditions. At present this discrep- 
ancy is unresolved. However it might be related to the in- 
herent assumption of a circular polar cap in equation 10. 
Indeed Table 2 does show different values of 6pc for the day- 
and nightside indicating non-circular polar caps. 

Finally in scenario N we show the reference scenario of 
the present Sun-Earth interaction. 

6. Discussion and conclusion 

We have carried out a study of the young Sun's solar wind 
and its interaction with Earth's paleomagnetosphere. Using 
numerical MHD simulations we have modeled the solar wind 
and evaluated its sensitivity to sunspot placement, strength 
of sunspots' magnetic fleld, strength of dipole fleld, coronal 
base density and rotation period. We have parameterized 
our solar wind model so as to agree with present solar wind 
data taken at 1 AU. We have also analyzed how sensitive the 
solar mass flux distribution is to the aforementioned param- 
eters. In comparing these findings to scaling relationships 
found in the literature we have determined that our mass 
flux estimates for the young Sun are approximately an order 
of magnitude less, with the exception of those solar wind so- 
lutions that have a higher coronal base density as their initial 
condition. We have found that the coronal base density is 
the controlling factor in mass flux magnitude, whereas mass 
flux distribution seems mainly affected by sunspot and back- 
ground dipole magnetic fleld strength. Sunspot latitudinal 
placement seems to affect mass flux slightly but leaves mass 
flux distribution unaffected. 

Given literature estimates as well as observations wc con- 
clude that simulations of the young Sun's solar wind likely 
require a high coronal base density. 

Based on our solar wind simulations we have established 
a number of solar wind scenarios with which to force the 
Earth's paleomagnetosphere. In these scenarios we varied 

parameters such as the Earth's dipole field strength and tilt, 
the solar wind speed, density, IMF orientation and strength. 
We have found that for a solar wind with high density, strong 
southward IMF and very weak dipole field the magnetopause 
stand-off distance can be as small as 4.25 Re- For a strong 
dipole field this increases to 8.5 Re and up to 9 Re for a 
northward IMF. 

It is clear from our simulations that the young Sun's solar 
wind will have had easier access to Earth's paleoatmosphere, 
especially in the case of a relatively weaker, with respect to 
the present, magnetic field: compare for instance scenarios 
L and N in Figure 7. As such, assuming a weak terrestrial 
magnetic field, our findings support those made by Tarduno 
et al. [2010], with the exception of their polar cap calcu- 
lation. Using an uncomplicated power-law based approach 
they found a greatly compressed paleomagnetosphere as well 
yielding obvious implications for heating and subsequent ex- 
pansion of Earth's exosphere with important ramifications 
for Earth's early atmospheric evolution. 
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Wc note however that convergence of the stellar mass loss 
rate scaling relationships by Wood et al. [2002, 2005], which 
Tarduno et al. [2010] based some of their conclusions on, 
with our solar wind simulations for high coronal base densi- 
ties is not conclusive evidence that either are correct. Mea- 
surements of stellar coronal base densities, for young solar 
analogs, are sparse and the myriad uncertainties that accom- 
pany the aforementioned scaling relationships are well doc- 
umented. Until more and better observations are made of 
stellar mass loss rates, stellar coronal base densities and/or 
the number of assumptions and uncertainties are reduced in 
the scaling relationships, our simulations of the paleomag- 
nctosphoric response to the young Sun's solar wind give the 
most realistic results to date. 

We emphasize that this study was carried out for the qui- 
escent Sun and modeled only steady-state interactions. As 
such, transient, dynamic interactions between the Sun and 
Earth's paleomagnetosphere have not been modeled. As 
noted by Tarduno et al. [2010], the magnetosphere could be 
more compressed during strong CME events, that are ex- 
pected to be much more frequent and energetic on a young 
solar-like star, than our steady-state scenarios suggest. In 
this context, it will be important to understand the relative 
contributions of steady solar winds and impulsive events to 
the mass loss experienced by more active solar analogs. On 
the other hand, assuming a close correlation between the 
CME source region and the distribution of active regions on 
the solar disk [Gopalswamy et al., 2008], wc cannot dismiss 
the possibility that most CMEs at that tirnc came from high 
latitudes on the disk and therefore the number of CMEs ac- 
tually reaching the Earth was in fact smaller than at the 
present time. 

Having obtained a range of magnetosphere sizes based on 
physically realistic solar wind scenarios, this study provides 
a well-grounded framework within which to undertake fur- 
ther paleomagnetospheric and paleoatmosphcric evolution 
studies. 
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Figure 1. Magnetograms used in Figure 2. MDI 
magnetograms used in sensitivity analysis of solar wind 
speed, density, magnetic field strength and mass flux to 
sunspot placement. Top-left: original magnetogram for 
a solar maximum (CR1958); top-right: magnetogram for 
CR1958 with spots placed at ±30° latitude; bottom-left: 
magnetogram for CR1958 with spots placed at ±60° lat- 
itude; bottom-right: original magnetogram for a solar 
minimum (CR2074). 
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Figure 2. Seven SC solar wind solutions at 1 AU 
along Earth's orbit with a sunspot magnetic field scal- 
ing factor of 2, a dipole magnetic field scaling factor of 2 
and for a 27 day rotation period. The period covers 28 
days. The panels show plasma speed u (top-left), mag- 
netic field strength B (top-right), plasma number den- 
sity n (bottom- left) and mass fiux (bottom-right) . In the 
legend 'shift' refers to which latitude the sunspots have 
been shifted to. The solution called 'up' refers to the so- 
lar wind solution evaluated 15'' above Earth's orbit and 
'down' below. The 'hd' solution refers to a coronal base 
density of 10^" cm~^ instead of the regular 5x10* cm~^. 
We compare the solar wind solutions to ACE data at 1 
AU for CR1958 (from January 1, 2000, to January 29, 
2000) as well as the average values of ACE data for the 
year 2000. 
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Figure 3. Mass flux maps at 1 AU where we vary the 
sunspot latitudinal placement. All solutions are 'sf2df2', 
for a 27 day rotation period and have {Mq yT~^ deg~^) as 
unit. Note that 'shiftGO' is the same scenario as 'sf2df2' 
in Figure 6. 
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Table 1. Average solar wind parameters for different so- 
lar wind solutions extracted along the trajectory of the ACE 
spacecraft, which is situated at the LI Earth-Sun Lagrange 
point, at ^ 1 AU. The mass flux decreases for stronger dipole 
magnetic field and increases for stronger sunspot magnetic 
field since the stronger spots close down open flux. Hence 
the fast wind is reduced and the corona is dominated by slow 
wind, which has a higher density, which seems to dominate 
the overall behavior of the mass flux. We did not calculate 
the mass flux for the "up" and "down" solutions, "hd" refers 
to a coronal base density of 10^" cm~^ instead of the regular 
5 X 10* cm~3. 
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Table 2. 


Input 


and output 


parameters 


for the 


solar 







wind/palcomagnctosphcrc scenarios used in the GM MHD 
simulations. Ouput parameters arc: Magnctopausc subsolar 
stand-off distance; Magnctopausc dawn-dusk flank distance; 
Magnctopausc north-south flank distance; Polar cap magnetic 
co-latitude (dayside / nightsidc) ; Polar cap area normalized 
by present Earth scenario N. 

°: North flank distance. South flank distance is 11.25 Re- 
North flank distance. South flank distance is 11 Rg. North 
flank distance. South flank distance is 14.25 Re- Present 
value of polar cap co-latitude: 18.1° [Tarduno et al., 2010] 
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Figure 4. Magnetograms used in Figure 5. Modified 
MDI magnetograms (CR1958) used as input to the Solar 
Corona module of SWMF to find a solution for the solar 
wind. We vary the sunspot magnetic field strength (sf) 
and the background dipole magnetic field strength (df ) by 
a scaling factor. All magnetograms have sunspots placed 
at ±60° latitude and are for a 27 day rotation. Top-left: 
SF2DF5; top-right: SF2DF10; bottom-left: SF5DF2; 
bottom-right: SF10DF2. 



X - 14 



STERENBORG ET AL: YOUNG SUN-EARTH INTERACTION 




"I . . . < . . 1 10""! • • • • ' • 1 

100 200 300 400 500 600 700 100 200 300 400 500 600 700 

time (hr) time (hr) 



Figure 5. Eight SC solar wind solutions at 1 AU where 
we vary the sunspot magnetic field scaling factor (sf ) and 
the dipole magnetic field scaling factor (df). All solutions 
are for a magnetogram with sunspots placed at 60 degrees 
latitude and for a 27 day rotation period unless otherwise 
indicated by '15d' meaning a 15 day rotation. As a ref- 
erence we show ACE data at 1 AU for CR1958 as well as 
the average values of ACE data for 2000. 
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Figure 6. Mass flux maps at 1 AU where we vary the 
sunspot magnetic field scaling factor 'sf , right column, 
and the dipole magnetic field scaling factor 'df, left col- 
umn. All solutions except for the last are for a 27 day 
rotation period unless otherwise indicated by '15d' mean- 
ing a 15 day rotation, and have {Mq yr~^ deg~^) as unit. 
Note that the last row has different limits. 
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Figure 7. GM solutions for various parameter choices. 
Shown are cross-sections of the magnetosphere along the 
noon-midnight plane. The relevant parameters are in- 
dicated above each figure: 'plasma spocd/plasma num- 
ber density/IMF z-componcnt/Earth's dipolc magnetic 
ficld/dipolc tilt. The top row varies Earth's magnetic 
field from 133% to 25% of the present dipolc. The sec- 
ond row varies dipole tilt from 0° to 90° ('pole-on'). 
The third row shows from left to right a solution based 
on 'sf2df2' - 'sf5df2' with a weak dipole field, a solu- 
tion with northward b^J^-^ , and a solution based on the 
'sf2dfl0_15d_hd' solar wind solution. The bottom row 
shows respectively solutions based on: 'sf2dfl0-15d_hd' 
with a stronger Earth magnetic field, 'sfl0df2_15d_hd' 
and a solution representing present Sun-Earth interac- 
tion. 



